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Abstract. We give a brief survey of recent results concerning almost diagonal- 
ization of pseudodifferential operators via Gabor frames. Moreover, we show new 
connections between symbols with Gevrey, analytic or ultra-analityc regularity and 
time-frequency analysis of the corresponding pseudodifferential operators. 



1. Introduction 

The objective of this paper is to report on recent progress concerning time-frequency 
analysis of pseudodifferential operators with smooth symbols. We discuss in particular 
connections between, on one side, symbols in the classical Gevrey or (ultra-)analytic 
spaces, and, on the other, (ultra-) modulation spaces, Gabor frames, which are the 
framework of time-frequency analysis. 

Modulation spaces having moderate weights were introduced by Feichtinger in 1983, 
and since then have been extensively studied by Feichtinger and many other authors, 
cf. in particular the textbook [T7], to which we address for further references. We 
recall their definition and main properties in the subsequent Section [2.31 but the main 
insight to understand their introduction being to consider the decay property of a 
function with respect to the space variable and the variable of its Fourier transform 
simultaneously. The time-frequency representation employed for their definition is the 
so-called short-time Fourier transform (STFT), whose building blocks are the linear 
operators of translation and modulation (so-called time-frequency shifts) given by 

T,fi-) = f{--x) and M5/(.) = e2-^«7(-), x, ^ e M'^. 

For z = (x,^) we shall also write iT{z)f = M^T^f- Indeed, if is a non-zero window 
function in the Schwartz class 5(M'^), then the short-time Fourier transform (STFT) 
of a a function/tempered distribution / with respect to the window g is given by 

(1) Vgf ix,0 = {f,M^T,g) = [ fit) git -x) e'^-?* . 

In short: we say that / is in the modulation space Mm'^iW^), 1 < p,q < oo, where m 
is a moderate weight function (hence with at most exponential growth at infinity) if 
Vgf G L^iM.'^'^). So it is clear that the bigger the growth of m the faster the decay of 
Vgf is at infinity. 
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Starting from signal analysis and quantum mechanics, these spaces have then been 
proved to be the right spaces for symbols of pseudo differential operators (see [I9l [25l 
EU [ini |32] and references therein), in particular localization operators (cf. e.g. [9]), 
Fourier multipliers ([HE^) recently Fourier integral operators of Schrodinger type 
O El El El [H], providing a good framework for the study of PDE's [T2] . 

If we limit the use of time-frequency shifts on a lattice A = aZ'^ X PZ"^ c M^"^ and 
fix a window function g G L^(]R'^), then the sequence Q{g,K) = {gm,n = MnTmg, 
(m, Ti) G A} forms a so-called Gabor system. The set G{g, A) is a Gabor frame, if there 
exist constants A, B > such that 

^Il/ll2< E K/.5™,n)P<i?||/||i, V/GL2(M'^). 

(r?i,n)£A 

For details see the next Section 2, where beside modulation spaces and Gabor frames, 
we treat classes of symbols and Gelfand-Shilov functions, which represent the natural 
window functions in our context. 

The first paper which uses Gabor frames to approximately diagonalize pseudodiffe- 
rential operators is [29]. Later, motivated by the study of Sjostrand [31j, Grochenig 
proved an almost diagonalization for symbols in the Sjostrand class, that is the mod- 
ulation space M°°'^, whose further generalization in |21) is recalled in Theorem 13.11 
below. What is remarkable in Theorem 13.11 is that the rate of decay in the almost di- 
agonalization, expressed in terms of Wiener amalgam spaces, characterizes the symbol 
class, which reveals to be a class of modulation spaces. The discrete representation of a 
pseudodifferential operator via Gabor matrix is one of the key ingredient for numerical 
applications, as shown in for the Gabor matrix of Schrodinger propagators. 

Motivated by these characterizations and by numerical applications, in the recent 
work [13] we study diagonalization of pseudodifferential operators with symbols of 
Gevrey, analytic and ultra-analytic type, and show that Gabor frames surprisingly 
reveal to be an efficient tool for representing solutions to hyperbolic and parabolic- 
type differential equations with constant coefficients. 

Here we are mainly focused on the theoretical aspects of the almost diagonaliza- 
tion: as a preliminary step we establish the equivalence between Gevrey, analytic and 
ultra-analytic regularity of symbols, and their membership to modulation spaces with 
sub-exponential, exponential and super-exponential weights, respectively. For these 
aspects, our results are strictly related to those of [3l]. The main result of the paper. 
Theorem 13.41 characterizes the corresponding classes of pseudodifferential operators 
in terms of exponential decay of the almost diagonalization. In the last Section 4, we 
report on the sparsity result obtained in [13], which motivates this study for numerical 
applications. 

Notations. The Schwartz class is denoted by 5(M'^), the space of tempered distribu- 
tions by S'{R'^). We use the brackets {f,g) to denote the extension to 5'(M'^) x 5(M'^) 
of the inner product {f,g) = J f{t)g(fjdt on L^(R'^). 

Euclidean norm of x E M'^ is given by |x| = (xf + • • • + x^)^^^ , and (x) = (1 + 
|3;|2ji/2_ \Ye write xy = x ■ y for the scalar product on M*^, for x,y € M*^. 



A FEW REMARKS ON TIME-FREQUENCY ANALYSIS 3 

For multiindices a, P £ N'^, we have \a\ = ai + ■ ■ ■ + ad, al = ail ■ ■ ■ a^!, = 

• • • <^ X G and, for /3 < a, i.e. /3, < a,, j G {1, 2, . . . , d}, (^) = {fj (^^^). 

The letter C denotes a positive constant, not necessarily the same at every appearance. 
The operator of partial differentiation d is given by 

oa _ oo _ gai _ _ _ Q«tj 
^ ^Xl x^ 

for all multiindices a £ N'^ and all x = (xi, . . . , x^) G R*^. If / and g are smooth enough, 
then the Leibnitz formula holds 



Observe 



I3<a 

and 

(3) |a|!<dl°la!. 

The Fourier transform is normalized to be f{^) = J-f{Cj = f /(t)e~^'^**^dt. 

For < p < oo and a weight m, the space £m(A) is the (quasi- )Banach space of 
sequences a = {ax}xeA on a lattice A, such that 

i/p 

\a\\£P 



:= |aAr"i(A)f^ < oo 

\AeA / 



(with obvious changes when p = oo). 

We shall use the notation A < B to express the inequality A < cB for a suitable 
constant c > 0, and A^ B for the equivalence c~^B < A < cB. 



2. TiME-FREQUENCY ANALYSIS AND GeLFAND-ShILOV SPACES 

2.1. Gelfand-Shilov Spaces. The Schwartz class 5(M'^) does not give enough infor- 
mation about how fast a function / G 5(M'^) and its derivatives decay at infinity. This 
is the main motivation to use subspaces of the Schwartz class, so-called Gelfand-Shilov 
type spaces, introduced in [TB] and now available also in a textbook [28j. Let us recall 
their definition and some of their properties. 

Definition 2.1. Let there be given s,r >0. A function f G 5(M'^) is in the Gelfand- 
Shilov type space S^{W^) if there exist constants A,B>0 such that 

(4) |x"9^/(x)| < Al"lsl^l(a!)'■(/3!)^ o./JgN"^. 

The space S^(K'^) is nontrivial if and only ifr + s>l, orr-|-s = l and r,s > jl6j . 

1 /2 r1 

So the smallest nontrivial space with r = s is provided by S-^^^^{M. ). Every function of 
the type P(x)e~"l^l^, with a > and P{x) polynomial on M*^, is in SljliR'^). Observe 
that (M^) C S^'^{R'^) for si < S2 and n < r2. Moreover, if / G S^{R'^), for every 
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5, 7 G N'^, x^d'^f G S^{R'^). The action of the Fourier transform on 5^(M'^) interchanges 
the indices s and r, as explained in the following theorem. 

Theorem 2.2. For f € S{R'^) we have f G S^{R'^) if and only iff G S^' 



Therefore for s = r the spaces S'KM ) are invariant under the action of the Fourier 
transform. 

Theorem 2.3. Assumes > 0, r > 0,s+r > 1. For f G tS(M'^), the following conditions 
are equivalent: 

a) / G Sm") . 

b) There exist constants A,B > 0, such that 

\\xy\\L^<A\"\{a\Y and < i^'^'W, a,f3eN^ 

c) There exist constants A,B > 0, such that 

\\xV\\L^<A\"^{a\Y and Wd" f < B^f'^ {/Siy , a,/3GN'^. 

d) There exist constants h,k > 0, such that 

||/e''''''^''''||L°° < oo and ||/e''l^l^''''||L°o < oo. 

A suitable window class for weighted modulation spaces (see Definition 12.81 below) 
is the Gelfand-Shilov space Sj;(M'^), consisting of functions / G 5(]R'^) such that for 
every constant A > and B > 

(5) \x''d^f{x)\<A\'^\B\^\a\/3\, a,^GN'^. 

We have S'KM'^) C Il}(]R'^) C 5j(M'^) for every s < 1. Observe that the characterization 
of Theorem 1 2 . 3 1 can be adapted to S]^(]R'^) by replacing the words "there exist" by "for 
every" and taking r = s = 1. Similarly, one can define Sg(R'^) for any r > 0,s > 0, 
with r + s > 1. 

Let us underline the following property, proved in [28, Proposition 6.1.5] and [THj 
Proposition 2.4], exhibiting two equivalent ways of expressing the exponential decay 
of a continuous function / on M'^. 

Proposition 2.4. Consider r > and let h be a continuous function on W^. Then 
the following conditions are equivalent: 

(i) There exists a constant e > such that 

(6) <e-^l^l^ xGR"^, 

(ii) There exists a constant C > such that 

(7) |x"/i(x)| < Cl°l(a!f , xGR*^, aGN'^. 

Remark 2.5. To he precise, the relation between the constants e and C is as follows. 

frdV 

Assuming then ^ is satisfied with C = ( — j . Viceversa, ([7]) implies ([6]) for 

any e < r{dC)^~ . Also, it follows from the proof that the constant implicit in the 
notation < in ([6]) depends only on the corresponding one in ([7]) and viceversa. 



A FEW REMARKS ON TIME-FREQUENCY ANALYSIS 



5 



The strong dual spaces of 5^(M°') and S];(M'^) are the so-called spaces of tempered 
ultra-distributions, denoted by (S'^)'(IR'^) and (Sj)'(M'^), respectively. Notice that they 
contain the space of tempered distributions 5'(M'^). Moreover the spaces S^{W^) are 
nuclear spaces [26] , This provides a kernel theorem for Gelfand-Shilov spaces, cf. [35] . 

Theorem 2.6. There exists an isomorphism between the space of linear continuous 
maps T from S^{M.'^) to {S^y{W^) and (S^)' (M.'^'^) , which associates to every T a kernel 
Kt G (5^)'(M2'^) such that 

{Tu, v) = {Kt, V u), Vu, v G 5^(M'^). 

Kt is called the kernel ofT. 

2.2. Gabor frames and time-frequency representations. We recall the basic 
concepts of time-frequency analysis and refer the reader to |17J for the full details. 
Consider a distribution / G 5'(M'^) and a Schwartz function g G 5(M'^) \ {0} (the 
so-called window). The short-time Fourier transform (STFT) of / with respect to g is 
defined in ([1]). The short-time Fourier transform is well-defined whenever the bracket 
(•, •) makes sense for dual pairs of function or (ultra-)distribution spaces, in particular 
for / G S'iW^) and g G S{R'^), f,g (£ L^{R'^), f G (S})'(R'^) and g G Sj(]R'^) or 
/ G (5^)'(M'^) and g G 5^(M'^). 

Another time-frequency representation we shall use is the (cross-)Wigner distribu- 
tion of f,g £ L2(M«'), defined as 

(8) W{f, g){x, = j^j[x+\^g{x- i)e-2-*« dt. 

If we set g{t) = g{—t), then the relation between cross- Wigner distribution and short- 
time Fourier transform is provided by 

(9) W^(/,<7)(x,e) = 2V-"«Fg/(2x,2e). 

For the discrete description of function spaces and operators we use Gabor frames. 
Let A = AZ'^'^ with A G GL{2d, M) be a lattice of the time-frequency plane. The set 
of time-frequency shifts Q{g,A) = {7r{X)g : A G A} for a non-zero g G L^(M'^) is called 
a Gabor system. The set Gig,-^) is a Gabor frame, if there exist constants A,B>0 
such that 

(10) AWfWl < ^ \{f,7,{X)g)\' < BWfWl V/ G L\R'). 

AeA 

If (fTO|) is satisfied, then there exists a dual window 7 G L^(M'^), such that ^7(7, A) is a 
frame, and every / G L'^iR'^) has the frame expansions 

/ = j;(/,vr(A)5)^(A)7 = ^(A)7)vr(A)g 
AeA AeA 

with unconditional convergence in L^(M'^). 

Eventually, we list some results about time-frequency analysis of Gelfand-Shilov 
functions, cf. [H [22l [33] : 

(11) /,5 G St(R^), s > 1/2 ^ Vgf G StiR^^), 
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Uge SI{R'^), s > 1/2, then 

(12) / G 5f (M^) ^ \Vg{f){z)\ < e-^l^l'^' for some e > 0. 

Since we have not found a precise reference in the hterature, we also recaU and 
present the proof of the following properties concerning the Wigner distribution. 

Proposition 2.7. We have: 

(13) f,ge SliR''), s > 1/2 ^ Wif,g) G SHR^"), 

(14) f,g e ^ W if, g)ei:\iR^''). 

Proof. We prove (I13p . formula (jl4p follows by similar techniques. The assumption 
g G SI{R'^) trivially implies g G SKM'^) and ([H]) proves that the short-time Fourier 
transform Vgg is in S'KM^'^). To prove that $ = W{g,g) G S'KR^'^), we use the 
equivalence a) 44> c) in Theorem 12.31 and the connection between Vgg and the Wigner 
distribution in ([9]). First, let us prove the first inequality in Theorem 12.31 c). For 
z = {x, ^) G M^'^, we have 

\z^W{g,g){z)\=2''\z^Vi^g{2z)\=2''-\''\\{2zrVi^g{2z)\ 

< 2'^-l"l||z"yg<7||Loo < 2'^-l"UI"l(a!)^ < (I) " (aiy. 

The second inequality in Theorem 12. 31 c) requires more computations but the techniques 
are the same as before. Using Leibniz formula for the x-derivatives and the ^-derivatives 
separately, 

d^d^,W{9,g){x,0 = (af (2'^e^--«y55(2x,20) 

= 2'^dl^ (^Q(47rie)^"^e^^^^^2lTl92(1^g5)(2a;,2e) 

= 2^ E (0 E ((4vr.0'^-)(4vr.x)^^e--« 



2l7l+|53l5|3 57(y.^)(2:,,20. 



E 2H+l'^3|(4^,)|/3-7l+l^ 



^l</3-7 

X e'^-^-''ix^2g4^ia;5^|3 57(y„^)(22;, 2^). 

Now, we have 

\e-'^-''x'^d'^dl{V0){2x, 2i)\ < Al"l+l^l((/3 -J -6^ + 62 + 63+ j)\r 

= Al"l+l^l((/3 + a-25l)!)^ 
where we used 82 + 63 = a — 5i. Since from ([2]) 

, < 2^l^-^^6i\ < 2l^-Tl((25i)!)i/2 < 2l''-^l((25l)!)^ 

(/3-7-5i)! 
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together with Y.^<p (^) = 21'^! and Y.&^+S2+&z=a 5^\5l\5s\ = obtain the desired 

estimate 

\dtd^W{g, g){x,0\< ^"+^1 ((« + ■ 



2.3. Modulation Spaces. Modulation spaces measure the decay of the STFT on the 
time- frequency (phase space) plane and were introduced by Feichtinger in the 80's [15], 
for weight of sub-exponential growth at infinity. The study of weights of exponential 
growth at infinity was developed in |14t I33j. 

Weight Functions. In the sequel v will always be a continuous, positive, even, 
submultiplicative function (submultiplicative weight), i.e., f(0) = 1, v{z) = v{—z), 
and v{zi + Z2) < v{zi)v{z2), for all z,zi,Z2 € M^'^. Submultiplicativity implies that 
v{z) is dominated by an exponential function, i.e. 

(15) 3C,k>0 such that I < v{z) < Ce''^'^ , zeR^"^. 
For instance, weights of the form 

(16) v{z) = e'\'\\l + \z\y log^(e + \z\) 

where a, r, s > 0, < 6 < 1, satisfy the above conditions. 

We denote by ^Ay{M.'^'^) the space of v-moderate weights on M^'^; these are positive 
and measurable functions m satisfying m[z -|- C) < Cv{z)m{Q for every z,(^ £ M^*^. 

Definition 2.8. Given g G S];(]R'^), a weight function m G A1„(M^°'), and 1 < p,q < 
00, the modulation space Mm'^{M.'^) consists of all tempered ultra- distributions f G 
(Si)'(M'^) such that Vgf G L^fiM.'^'^) (weighted mixed -norm spaces). The norm on 



IS 



(17) II/IIm- = II^./IIl- = I / {/ \Vgf{x,CWm{x,0''dx\ d^ 




(obvious changes if p = 00 or q = 00). 

For f,g £ the above integral is convergent and thus S|(M'^) C M^'^(M'^), 

1 ^ P,q ^ 00, cf. with dense inclusion when p,q < 00, cf. [3j. When p = q, we 
simply write M^(M'^) instead of M^^(M'^). The spaces M^''(]R"') are Banach spaces 
and every nonzero g G M^(R'^) yields an equivalent norm in ()17p and so Mm'^{W^) is 
independent on the choice of g G M^(M'^). 

We observe that these properties of modulation spaces do depend on the fact that 
the weight functions involved have at most exponential growth at infinity. Indeed, as 
well known, if we consider super exponential weights at infinity, say, e.g. m{z) = e'^'^l , 
with A; > and 1 < 6 < 2, then the related modulation spaces Mm'^{M.'^) can still be 
defined by taking Gelfand-Shilov windows, but their definition depends on the choice 
of the window g, because of the loss of sub-multiplicativity of the weight. 
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2.4. Ultra-Modulation Spaces. In this section we present a new definition of mod- 
ulation spaces, whicli considers also weights of super-exponential growth at infinity (of 
a very particular form). For related constructions, we refer to [3l|22] and more recently 
to |34j- Let us first introduce the weights: 



(18) Ws,e{z) := e^\'\\ z € M"', s > 1/2, e > 0, 

and the weight class AA(M2'^) = {ms^e{x,C) = (1 (8) Ws,e){x,^) = Ws,e{C), x,^, gW^, e > 
0, s > 1/2}, which we may identify with the set (0, +oo) x [1/2, -|-oo) of the correspond- 
ing e, s parameters. Thus, a weight mg^e £ J\f{M?'^) grows faster than exponentially at 
infinity in the frequency variable whenever 1/2 < s < 1. Notice that the bound 
s = 1/2 is admitted. We use the class AA(M^'^) as weight class for modulation spaces. 
We limit our study to weights in the frequency variable for simplicity. On the other 
hand, observe that a limit in enlarging the definition to the space variables is imposed 
by Hardy's theorem: if m(z) > Ce'^'^l , for z = (x,^) G M^*^ and some c > vr/2, then 
the corresponding modulation spaces are trivial [22j . 

Definition 2.9. Let e > 0, s > 1/2 and rus^e G M{R'^'^). Consider a non-zero window 
function g in Sg{W^). For 1 < p,q < oo, we define Mm^.^^gi'^'^) the subspace of 
f € (S];)'(M'^) such that the integrals in (fT7|) are finite (with obvious changes if either 
p = oo or q = oo). 

Notice that: 

(i) If s > 1/2, /, 5 G Sy'^{W^), the integral in (fT7|) is convergent thanks to (fTT]l and 
Theorem 12.31 item d). Indeed, we can find h> such that || V^/e^'''^ ||loo < oo and 

\ <i/p 

\Vgfix,0\''ms4x,0''dx] di 

<C\\{VJ)e''\-\'\\L^ [ ([ \msAx,CWe-^P\^^'^^\"dxY%<oo. 

(ii) For nis^e & Mv(R'^'^), hence s > 1, then M^f^,g(M'^) is the subspace of ultra- 
distribution (Sj)'(R'^) defined in Definition 12.81 So we come back to the classical 
modulation spaces. This justify the same notation. 

(iii) The definition of Mms,e,g{R'^) may depend on the choice of the window function 
g. However, if s > 1, the definition of Mms,c,g{R'^) does not depend on g: the class of 
admissible windows can be enlarged to S'^'^) [S Proposition 1], so we recapture 
the standard definition [ 



2.5. Time-frequency analysis of Gevrey-analytic and ultra-analytic symbols. 

Here we present the results obtained in \13\ Section 3] , where the smoothness and the 
growth of a function / on M.'^ is characterized in terms of the decay of its STFT Vgf, 
for a suitable window g. 

Theorem 2.10. Consider s > 0, m G 7W^(M''), g G M^^^iR"^) \ {0} such that there 
exists Cg > 0, 

(19) l|9°<7||il(M.)<Ci"l(a!)^ a en''. 
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For f € C°°(M'^) the following conditions are equivalent: 
(i) There exists a constant Cf > such that 

(20) |9°/(x)| <m(a;)C7|"l(a!)^ x G M^, a G N^. 
(a) There exists a constant Cj^g > such that 

(21) mfix,^)\<mix)c\'^liair, (x, G a G N^. 
(Hi) There exists a constant e > such that 

(22) \Vgfix,0\<Hx)e-'^^^\ (x,e)GM2^aGN'^. 

We say that the function / is Gevrey if s > 1, analytic if s = 1, and ultra-analytic 
when s < 1. 

Remark 2.11. We observe that the assumption m G Aiv{^'^) in the previous theorem 
is essential to obtain the equivalence: indeed in the proof we exploit the v -moderateness 
of m. 

For simphcity, from now on we shah assume m = v = 1. 

A natural question is whether we may find window functions satisfying (jl9p . To 
find the answer, we recall the following characterization of Gelfand-Shilov spaces. 

Proposition 2.12. Let g G S{W^). We have g G S^^{W^), with s,r > 0, r + s > 1, if 

and only if there exist constants ^ > 0, e > such that 

|a"5(x)| < ^l"l(a!)"e-^l^l\ X G M"^, a G n'^. 

We have g G S];(M'^) if and only if, for every ^ > 0, e > 0, 

|a"5(x)| < Al"la!e-^l^'l, X G M"^, a G N'^. 

Hence every g G S'^(IR'^) with s>0, 0<r< l,s + r> 1, satisfies (1190 for every 
submultiplicative weight v (see ()15p ). The same holds true if 5 G and s > 1. 

3. Almost Diagonalization for Pseudodifferential operators 

In this section we first present the almost diagonalization for pseudodifferential 
operators having (ultra-) analytic symbols obtained in \12>\ Section 4]. This result can 
be seen as an extension of the almost diagonalization for pseudodifferential operators 
obtained in [21], where only the Gevrey-analytic case was discussed. 

Secondly, we exhibit a new characterization involving ultra-modulation spaces. 

Recall the Weyl form a'^ of a pseudodifferential operator (so-called Weyl operator 
or Weyl transform) with symbol (t(x,^) on M^*^, formally defined by 

(23) a-/(x) = j^^a ^) ^2.ii.-y), ^ ^y^^_ 

Using the Kernel Theorem for Gelfand-Shilov spaces (Theorem 12. 6p , we have a char- 
acterization of linear continuous operators T : S^{M.'^) — > {Sf.y{M.'^). In particular, any 
such operator T can be represented as a pseudodifferential operator in the Weyl form, 
with a G (5J)'(M^^). Thus we shall exhibit our results for Weyl operators. 
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It was proved in [181 [2T] that Gabor frames allow to discretize a continuous operator 
from iS(M'^) to S'{W^) into an infinite matrix that captures the properties of the orig- 
inal operator. In particular, the authors consider symbols in the modulation spaces 
^i^y^j-ii where is a continuous and submultiplicative weight function on M?'^ and 
j is the rotation on M^'^: 

j{zi, Z2) = {Z2, -Zl), (Zi, Z2) G K^'^. 

Fix a Banach algebra of sequences ^?(A), where A is a lattice on M^"^ with rela- 
tively compact fundamental domain Q containing the origin. Then a function H G 
L^^^(M^'^) belongs to the Wiener amalgam space W{£t)(M.'^'^) if the sequence h{X) := 
esssnpy^(z\^QH{u) is in il{A). The almost diagonalization for pseudo differential oper- 
ators in [18\ [2T] can be extended easily to modulation spaces with exponential weights 
v{z) = 6*^1^1 [T3]- Indeed one can use the same pattern of [2I5 Theorem 4.1], pro- 
vided that the symbol space of distributions 5'(M^'^) is replaced by the space of ultra- 
distributions (S})'(M^'^). Thus we rephrase |2H Theorem 4.1] in our context as follows. 

Theorem 3.1 ([13 [21]). Assume that Q{g,A) is a frame for L'^ {W^) with g G M^i(M^). 
Then the following statements are equivalent for a G (S^)'(M^'^).- 

(ii) There exists a function H G W{£f,){M.'^'^) such that 

(24) \{a^TT{z)g,7T{w)g)\<H{w-z), yw,zeR^''. 
(iii) There exists a sequence h G ^?(A) such that 

(25) \{a^7r{X)g,7r{fi)g)\<h{f,-X), VA,^G A. 

The Gabor kernel of a'": {a''"TT{z)g,7r{w)g) is referred as continuous Gabor matrix, 
whereas {a'^iT{\)g,iT{^)g) is the Gabor matrix of cr"". 

Our concern is in the case of a continuous Gabor matrix dominated by a function 
H G W{l'^){E?'^). In this special case, Theorem 13.11 can be simplified as follows. 

Theorem 3.2. Assume that Q{g,K) is a frame for L'^{W^) with g G Ml{W^). Then 
the following statements are equivalent for a G (S];)'(M^'^); 

r^;^eM-„„^._,(M2rf). 

(ii)The continuous Gabor matrix satisfies 

(26) \{a^TT{z)g,7T{w)g)\<v{w-z), yw,zeR^''. 
(iii) The Gabor matrix satisfies 

(27) \{a'"7r{X)g,7r{fx)g)\<v{fi-X), VA,/zGA. 

Our main aim here is to give a counterpart of Theorem 13.21 for weights of super- 
exponential growth. First, let us survey the new results of [13^ Section 4] which provide 
the arguments for the characterization of Theorem 13.41 below. 

The crucial relation between the action of the Weyl operator a'^ on time-frequency 
shifts and the short-time Fourier transform of its symbol, contained in [18[ Lemma 3.1] 
can now be extended to Gelfand-Shilov spaces and their dual spaces as follows. 



Lemma 3.3. Consider s > 1/2, g E S: 

(28) \{a^7r{z)g,7r{w)g)\ 
and 

(29) \VMu,v)\ 
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,^ = W{g,g). Then, for a e {S! 

= \ Vq,a{u,v)\ 



. z + w 
l^$cr I — - — ,j{w - Z 



[a TT 



Proof. Since $ = W{g, g) G S|(IR2<i) for g € 5f (M'^) by the duality {a, 7t{u, v)^)(^siyxSi 
is well-defined so that the short-time Fourier transform V$cj(u, v) makes sense. The 
rest of the proof is analogous to [lU Lemma 3.1]. □ 

Given a pseudodifferential operator a"^, with smooth symbol a € C°°{M?'^), we 
exhibit the connection between the Gevrey, analytic or ultra-analytic regularity of a 
(as considered in Theorem 12.101 the dimension being now 2d), the (ultra-)modulation 
space a belongs to (cf. Definition 12. 9p . and the decay of the continuous Gabor matrix 
of a"". 

Recall the weights Ws,e defined in (jlSp and here used as functions over M?'^. 



Theorem 3.4. Let s > 1/2 and consider a window function g ^ Sg 
W{g,g). Then the following properties are equivalent for a G C°°(M^'^).' 

(i) There exists e > such that a £ M^^^ ^^^{R'^'^). 

(ii) There exists C > such that the symbol a satisfies 

(30) |5"cj(z)| < Cl"l(a!)^ G M^"', Va G N^'^. 
(Hi) There exists e > such that 

(31) \{a'"TT{z)g,-K{w)g)\ <Ws^^^{w - z), y z,w e 
Proof, (i) =^ {Hi). Assuming a G M^^^ ^ $(M^°'), that is 



Set $ 



r,2d 



sup |y$a(u,z;)|e^l^l' 



< oo 



and using (j29j) . we obtain the claim: 

\{a^7T{z)g,n{w)g)\ 



< sup \VMuJ{w-z))\<e-'\^-'\\ 



(Hi) =^ (i). Relation (l29]l and the decay assumption (fSTl) give 



\V^a{u, v) 



j (v)] g,TT[u + -j (v) g) 



1 

-e\v\ s 



(32) 
(33) 

The equivalences (ii) <^ (Hi) are proved in Theorem \13\ Theorem 4.2]. For sake of 
clarity, we give a brief sketch of this result. The techniques are similar to [18_l Theorem 
3.2]. 
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{ii) =^ {in). The window ^ = W{g,g) G SI{R'^'^), for g G 5|(M^) by and 
satisfies the assumptions of Theorem I2.10[ Hence, using the equivalence (f20]) 44> (f22]l . 
the assumption ()30p is equivalent to the following decay estimate of the corresponding 
short-time Fourier transform 

\Vi,a{u,v)\<e-^\''\\ u,v£R^'^, 

for a suitable e > 0, hence 

which combined with (|28p yields (iii). 

{Hi) (n). We use relation (j29p and the decay assumption (j3ip again, which give 
(j32p and using the equivalence (j20p <;=^ (j22p we obtain the claim. □ 

Remark 3.5. a) We observe that the constant e > 0, which depends on the choice 
of g ^ Sg{W^), is the same in {i) and {Hi). Whereas the link between e > and the 
constant C > in {ii) is specified in Remark \2.5[ 

b) If we consider s > 1, that is the symbol a is a Gevrey or an analytic symbol, then 
the equivalence {i) 44> {ii) follows by combining Theorems \ 3.^ and [iS', Theorem 4.2]. 

Of course the estimate (f3T]l implies the discrete analog (|34|) below. The vice versa is 
not obvious and requires the existence of a Gabor frame Q{g,A) having g G Sy'^{W^) 

and a dual window 7 G Sy2{^ ) well. We call such a Gabor frame a Gabor super- 
frame. The existence of Gabor super-frames is due to a result obtained by Grochenig 
and Lyubarskii in |20) . They find sufficient conditions on the lattice A = AZ"^, A G 
GL(2,M), such that g = Yl^=o'^kHk, with Hi^ Hermite function, forms a Gabor frame 
G{g,A). Besides they prove the existence of dual windows 7 that belong to the space 
S^^^C^) (cf- |20l Lemma 4.4]). This theory transfers to the d-dimensional case simply 

by taking tensor products g = gi®- ■ ■'S>gd G {^ ) of windows as above, which define 
a Gabor frame on the lattice Ai x • • • x and possess a dual window 7 = 71 18) • • • (X) 7^ 
in the same space G S'^^gC^^'^)- 

The Gabor super- frames allow the discretization of the kernel in ()3ip . 

Theorem 3.6. Let Q{g,A) a Gabor super-frame for Gonsider s > 1/2 and a 

symbol a G C°°(IR^'^). Then the following properties are equivalent: 

(i) There exists e > such that the estimate (j3ip holds. 

(ii) There exists e > such that 

(34) |(^7'"7r(^)g,7r(A)g)| <u;,,_e(A-/i), VA,^gA. 

Proof. Let us sketch {ii) =^ {i). The pattern of [181 Theorem 3.2] can be adapted to 

1 /2 r7 

this proof by using a Gabor super-frame Q{g,A), with a dual window 7 G Sy2{^ ) 

and the following property of the weights (fT8]l { fl3\ Lemma 3.5]): 

(35) 

r w,^ ^ , . ( Ws-e{X), for s>l 

{ws,-.*Ws,-.){X):=}_^^WsM>^-u)wsM^)<^ ^., _^2_./.(A), for i < . < 1. 
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Let Q be a symmetric relatively compact fundamental domain of the lattice A C 
Given z G M^'^, we can write them uniquely as w = X + u, z = + u' , for \, fx £ A 
and u, u' S Q. Using the Gabor reproducing formula for the time-frequency shift 

1 /2 r] 

TT{u)g G Sy^iW^) we can write 

'^{u)9 = ^{j^{u)9,'^iy)l)T^{v)9- 

Inserting the prior expansions in the assumption (j34p . 
|(a"'7r(;U + n')5,vr(A + u)5)| 

< K<T"^(/" + ^05,vr(A + i/)g)||Mn05,4^')7)IIMn)5,vr(z.)7)| 
i/.i/'eA 

< E -^i ^^^X"^"' ) e-l^+-^-'l- \V,9{.' - u')\\V,g{u - u)\. 

Since the window functions 5,7 are both in Sy'^{M.'^), the STFT Vyg is in Sy"^ 
due to ifTT]) . Thus there exists h > such that |V^(7(z)| < e"'^'^'^, for every 2; G 



p2d 



Inserting this estimate in the previous majorization and using ()35p repeatedly, 
|(a-7r(;U + n')5,vr(A + ^x)5)|< ^ ^~ 



< ^ ^-b{\\ + l^-fj.-u'\s +\u'\s) 

= Ws,b * WsM * 'Ws,b)(A - fx) 



< e 



-e\X—^j,\ ~' 



for a suitable e > 0. 

If 2; G M^'^ and w = X+u, z = /i+u', A, ;U G A, u, u' G Q, then A—// = w — z+u' —u 
and u' — u Q — Q, which is a relatively compact set, thus 



(36) e-^'l^-'^l^ < sup e-^>-"+"l^ < e-^>-"l' 

ueQ-Q 

This gives the desired implication. 



4. Sparsity of the Gabor matrix 

The operators in Theorem 13.41 enjoy a fundamental sparsity property. Indeed, 
let Q{g,A) be a Gabor super- frame for L'^{R'^). Then, as we saw, 

(37) \{a^7r{fi)g,7riX)g)\<Ce-'\^~^'\\ VA,/iGA, 

with suitable constants C > 0, e > 0. This gives at once an exponential-type sparsity, 
which is proved in [13, Proposition 4.5] (we refer to [21 [23] for the more standard notion 
of super-polynomial sparsity) and recalled below. 
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Proposition 4.1. Let the Gabor matrix {a^TT(fi)g,iT{X)g) satisfy ()37p . Then it is 
sparse in the following sense. Let a be any column or raw of the matrix, and let \a\n 
be the n-largest entry of the sequence a. Then, \a\n satisfies 



\a\n < Ce-^"''= , n G N 
for some constants C > 0, e > 0. 
Proof. By a discrete analog of Proposition 12.41 it suffices to prove that 

n"|a|n < C"+^(Q!)2'^^ Q € N. 

On the other hand we have 

"n^ ■ \a\n < 

for every < p < oo. Hence by (j37p and setting p = 1/a we obtain 

n"|a|„ < {j^e-'^^^-^^'^y = (^e-^fl^l^)^ 
aga aga 

Let Q be a fundamental domain of the lattice A. Then if x G A + Q, A G A, we have 
\x\ < |A| + Co, therefore < Ci(|A|i/" + 1). Hence 

E ,i [• ,1 [• r+oo 1 

g-.p|A|. <^ / e'^p\A^dx= da e-'PP^ p'^'^-^dp 

Css ^_t^2ds-,^^ _ C3sT{2ds) _ C, 



Finally, by Stirling's formula, 

p p 
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